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RESEARCH MEMORANDUM 


TIME-VECTOR DETERMINED LATERAL DERIVATIVES OF A 
SWEPT-WING FIGHTER-TYPE AIRPLANE WITH THREE 
DIFFERENT VERTICAL TAIIR AT MACH 
NUMBERS BETWEEN O. 7 O AND 1.1+8 
By Chester H. Wolowicz 


SUMMARY 


As part of the flight research program conducted on a swept-wing 
fighter-type airplane, rudder-pulse maneuvers were performed at altitudes 
from 30,000 to 1+3,000 feet over a Mach number range of O. 7 I to 1.1+8 to 
determine the lateral stability characteristics relative to the stability 
axes, in general, and the lateral derivative characteristics, in partic- 
ular. The time-vector method of analysis was used. Four configurations 
were employed in the investigation. Three configurations involved three 
different vertical tails with varying aspect ratio or area, or both. The 
fourth configuration employed a large tail, which had been used in the 
third configuration, and an extension of the wing tips. 


The time-vector method of analysis is capable of producing good 
values of the lateral derivatives Cy^, C n ^, C 1 , and Cj providing 

the damping ratio is less than approximately 0.3* Reliable values of 
lateral derivatives (Cn r - C n p) are difficult to determine because of the 

sensitivity of this quantity to other factors. The expected effects of 
increasing vertical-tail size, resulting in increased magnitudes of C n p, 

Cm, and , were realized. The addition of wing-tip extensions had 

small effects, except for a fairly large increase in the magnitude of the 
damping-in-roll derivative Cm. Theoretically calculated derivatives 

showed fair to good agreement with flight results in the subsonic range 
with the exception of high angle-of -attack values of (Cn r - Cnp) deriv- 
atives. Wind-tunnel data for the static derivatives for a Mach number 
of 1 . 1 + 1 , when corrected for torsional flexibility and air-intake effects 
of the jet engine, showed good agreement with flight results. The experi- 


mental rate of decrease in the magnitudes of C : 


v 


C 7 , and 

Zp> 


C 7 with 

b-n 


Mach number at Mach numbers greater than 1.25 was larger than estimated. 
This increased rate of decrease in magnitudes appears to be the result 
of possible shock wave and flow interference at the wing tips. 
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INTRODUCTION 


This paper presents an application of the time-vector method of 
analysis in the determination of static and dynamic lateral derivatives 
of a 45° swept -wing fighter -type airplane. Details of the application of 
the time-vector method of analysis employed are also included in this 
paper as are some considerations of the limitations of the method. All 
data were obtained as part of a comprehensive investigation, conducted 
at the NACA High-Speed Flight Station at Edwards , Calif., of the lateral 
characteristics of this airplane. 

The quantity of data obtained from the flight test program provided 
the first opportunity to perform a fairly detailed investigation of the 
stability characteristics of an airplane in the transonic and supersonic 
regions and to provide some comparison with available wind-tunnel data 
(ref. l) . Previous reports have presented the directional stability as 
determined by simple relationships (ref. 2) and the results of* roll cou- 
pling investigations (refs. 3 to 6) . 

The flight rudder-pulse data for the determination of the lateral 
stability characteristics were obtained for four conf igurations . Three 
configurations employed the original wing and three different vertical- 
tail areas (original, extended, and large) while the fourth configuration 
employed an extended wing and the large tail. The data were obtained at 
altitudes between 30,000 and 43 >000 feet over a Mach number range extending 
to 1.48. Most of the tests were performed at a nominal value of 1 g load 
factor, but for a few tests at nominal Mach numbers of 0.83 and 1.14 load 
factors within the range of 0.5g to 1.8g were used to investigate the 
influence of angle of attack on the lateral stability characteristics. 

The results of the analysis of the data are compared with available 
wind-tunnel data and calculated derivatives . 


SYMBOLS AND COEFFICIENTS 


The results of this investigation are referred to the stability 
system of axes, which is defined as an orthogonal system of axes inter- 
secting at the airplane center of gravity in which the Z-axis lies in the 
plane of symmetry and is perpendicular to the X-axis. The X-axis is in 
the plane of symmetry and is the projection in the XY-plane of the rela- 
tive airstream onto the XZ-plane of symmetry. The Y-axis is perpendicular 
to the plane of symmetry. 
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span. 

a n 

at 

at x 

b 

°io 

Cz 

% 


The coefficients are referred to the original wing area and wing 

normal acceleration, g units 
corrected transverse acceleration, g units 

indicated transverse acceleration uncorrected for instrument 
position, g units 

wing span, ft 

trim lg lift coefficient, W/qS 




C 7 • 


C 7 

£5: 




jn P 


rolling-moment coefficient, g^ngmoment 

qoD 

ac 7 

damping -in-roll derivative, — £, per radian 

c)£— 

2V 

rate of change of rolling-moment coefficient with yawing 

dC 7 

angular velocity factor, — p per radian 

% 


SC 


effective dihedral derivative, — -, per radian 

rate of change of rolling -moment coefficient with rate of 

Sc 

change of angle-of -sides lip factor, — per radian 

°2V 

rate of change of rolling-moment coefficient with respect 

Sc 

to control-surface displacement, per deg 

c6 

yawing -moment coefficient, Yswln g moment 

qSb 

rate of change of yawing-moment coefficient with rolling 

Sc 

angular velocity factor, — 2, per radian 
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rate of change of yawing -moment coefficient with yawing angular 

Sc 

velocity factor. — £, per radian 

°2V 

SC 

directional stability derivative , per radian 

dp 


rate of change of yawing-moment coefficient with rate of change 

Sc 

of angle-of-sideslip factor, — per radian 

dH 

°2V 


rate of change of yawing-moment coefficient with respect to 

dC n 

control-surface displacement, , per deg 


lateral-force coefficient. 


Lateral force 
qS 


rate of change of lateral-force coefficient with rolling 


dC 


angular velocity factor, — — , per radian 

^b 

°2V 


rate of change of lateral- force coefficient with yawing angular 

SC Y 

velocity factor, — — , per radian 

^rb 

2V 

dCy 

lateral-force derivative, — — , per radian 

dp 

rate of change of lateral-force coefficient with rate of change 

dCy 

of angle -of-s ides lip factor, — ri, per radian 

rate of change of lateral-force coefficient with respect to 

dC Y 

control-surface displacement, - — , per deg 

d& 


chord, ft 

mean aerodynamic chord, ft 
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g acceleration of gravity , ft/sec 2 

hp pressure altitude, ft 

lx moment of inertia of airplane about stability X-axis, 

IXq cos ^ t 1 + I Z q s in 2 il, slug-ft 2 

1x2, product of inertia referred to stability X- and Z-axes, 

_1 / 2 ( X Z 0 " %o) sin 2 t 1* slug-ft 2 

I^ moment of inertia of airplane about stability Z-axis, 

I Z 0 cos 2 i] + I X sin 2 r), slug-ft 2 

IXq;Iyq ; ^Zo momen ^ s °? inertia of air*plane about principal longitudinal, 
lateral, and vertical axes, respectively, slug-ft 2 

i^ incidence angle of horizontal tail, positive leading edge up, deg 

M Mach number 

m mass of airplane, W/g, slugs 

m a mass rate of air intake of jet engine, W a /g, slugs/sec 

P period of damped natural frequency of airplane, sec 

p T rolling angular velocity factor, pb/2V, radians 

q dynamic pressure, ^pV 2 , lb/sq ft 

r,p,p rate of change with time of \|r, cp, and p, respectively, 

radians/sec 

r,p rate of change with time of r and p, respectively, 

radians /sec 2 

r ! yawing angular velocity factor, rb/2V, radians 

S wing area, sq ft 

T-i /p time required for absolute value of transient oscillation to 
damp to half amplitude, sec 

t time, sec 


CONFIDENTIAL 


CONFIDENTIAL 


NACA EM H56C20 


V 

v 

W 

W a 

*a 

Xa t 

X P 

Za t 


a 

3 

Pi 


P' 

(*S) 


a 


airspeed, ft/sec 

transverse acceleration, ft/sec^ 

weight of airplane, lb 

weight rate of air intake of jet engine, lb/sec 

distance from center of gravity of airplane to air intake of 
jet engine (measured parallel to body X-axis), 25 . 1 ft 

distance from center of gravity to transverse accelerometer 
(measured parallel to body X-axis), positive when forward of 
center of gravity, 5-37 ft 

distance from center of gravity to sideslip vane (measured 
parallel to body X-axis), positive when forward of center of 
gravity, 30 ft 

distance from center of gravity to transverse accelerometer 
(measured perpendicular to body X-axis), positive when below 
center of gravity, -3-6 ft 

distance from center of gravity to sideslip vane (measured 
perpendicular to body X-axis), positive when below center 
of gravity, 2-5 ft 

angle of attack of airplane, angle between reference body 
X-axis and stability X-axis, deg 

corrected angle of sideslip, deg or radians 

indicated angle of sideslip, measured from relative airstream 
to X-axis, positive when X-axis is left of airstream, deg 

rate of change of sideslip factor, pb/2V, radians 

contribution of intake air of jet engine to directional 

-m a Vx a 


stability derivative. 


qSb 


-, per radian 


P/a 


contribution of intake air of jet engine to lateral-force 
-m a V 


derivative. 


qS 


-, per radian 


>a t 


total aileron deflection, positive when left aileron is down, 
deg 
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S r rudder deflection, positive when rudder deflected to left, deg 

e angle between reference body X-axis and principal X-axis, 

positive when reference axis is above principal axis at the 
nose, deg 

£ ratio of actual damping to critical damping 


n 


P 

T 

a 

a P 


°P 


$ 

*d 


angle of inclination of principal X-axis of airplane relative 
to stability X-axis, positive when principal X-axis is above 
stability axis at the nose, a - e, deg 

mass density of air, slugs/cu ft 

time parameter, m/pVS, sec 

angle of sidewash, radians 

rate of change of angle of sidewash with angle of sideslip, 

cto/dp 

rate of change of angle of sidewash with rolling angular 
velocity factor, 

^b 

2V 

phase angle, deg 
damping angle, deg 


cp angle of roll, positive when right wing moves down, radians 

\|r angle of yaw, positive when airplane turns to right, radians 

undamped natural frequency, radians/sec 




damped natural frequency. 




radians /sec 


(%v) F ’(%v) F ’ et0 


contribution of flexible, vertical tail to the 
lateral-force, directional-stability deriv- 
atives , etc . , respectively 



contribution of the rigid, vertical tail to the 
lateral-force, directional-stability deriv- 
atives, etc., respectively 
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change in the contribution of vertical tail to 
lateral-force , directional-stabi lity deriv- 
atives caused by flexibility of the vertical 
tail, etc., respectively 


The symbol | j | represents the absolute magnitude of a j quantity 
and is positive. When employed in an equation, the equation is consid- 
ered to be a vector equation. 


The phase angle of a vector j relative to another vector k is 
indicated by the .subscript The second subscript k is used as 

the reference. For example, in the expression = -150° the roll 

displacement vector lags the yaw displacement vector by I 5 O 0 . 


AIRPLANE 


The airplane is a fighter-type with a single turbojet engine equipped 
with an afterburner, a moderately low swept wing, and a low horizontal 
tail. A three-view drawing of the airplane with the original vertical 
tail, tail A, is shown in figure 1. Figure 1 also indicates the extended 
wing. A photograph of the airplane is shown in figure 2. 

The tests covered the following four configurations: 


Configuration 

Vertical tail 

Wing 

A 

Small (A) 

Original 

B 

Extended (b) 

Original 

C 

Large (c) 

Original 

D 

Large (C) 

Extended 


Figure 3 presents a photograph of tails A and C. Drawings of the three 
vertical tails are shown in figure 4. The same rudder was used on all 
tai Is . 


The airplane is equipped with automatic leading-edge slats in five 
interconnected segments. At subsonic speeds the slats generally started 
to open at 3° to 6 °. At supersonic speeds the slats generally remained 
closed at Mach numbers above 1.25 tor the angle-of-attack range of the 
tests . 


The physical characteristics of the various configurations are pre- 
sented in table I. The estimated variation with airplane weight of the 
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principal moments of inertia and inclination of the principal axes 
(fig- 5) is based on the manufacturer's estimate (ref. 7) for design 
weight and empty weight conditions. 


INSTRUMENTATION AND INSTRUMENT ACCURACY 


Standard NACA instruments were used to record airspeed, altitude, 
rolling and yawing velocities and accelerations, normal acceleration, 
transverse acceleration, angles of attack and sideslip, and rudder, 
aileron, and stabilizer positions. The airspeed, altitude, and angles 
of attack and sideslip were sensed on the nose boom. All records were 
synchronized at 0.1-second intervals by a common timing circuit. 

The turnmeters used to measure the angular velocities and acceler- 
ations were referenced to the body system of axes of the airplane and 
are considered accurate to within ±1.0 percent of scale range. Mounting 
direction errors were O.5 0 or less. 

The indicated normal and transverse accelerometer readings were 
corrected to the center of gravity. The accelerometers are considered 
accurate to within ±1.0 percent of scale range. 

Indicated sideslip angles and angles of attack, measured by vane- 
type pickups, were corrected for roll and yaw rate, and pitch-rate 
effects, respectively. The pickups were mass damped and had dynamically 
flat frequency-response characteristics over the frequency range of the 
airplane. The pickups are statically accurate to ±0.1°. 

The ranges, dynamic characteristics, and scale of recorded data for 
the angle of attack, sideslip, velocity, and acceleration instruments are: 


Function 

Range 

Scale of recorded 
data (per inch 
deflection) 

Undamped natural 
frequency, cps 

Damping 

ratio 

a, deg 

-20 to 40 

10.55 

8 

0.70 

?>, deg 

±32 

10.75 

8 

0.70 

r, radians/sec . . 

±0.5 

0.543 

10 to 12 

O.65 

r, radians/sec^ . . 

+1 

1.01 

8 

O.65 

p, radians/sec . . 

±4 

4.19 

20 

0.64 

p, radians/sec^ . . 

±7 

6.33 

7 

O.65 

a n , S units .... 

-1 to 7 

4.92 

a 19 

a 0.43 

a-j-, g units .... 

±1 

2.30 

b 13 

b 0.38 


a 3 0,000 feet 
Ho, 000 feet 
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Rudder, aileron, and stabilizer positions were measured by standard 
control-position transmitters linked directly to the control surfaces. 

The transmitter-recorder system had a flat dynamic response over the 
frequency range of the control movements encountered. The transmitters 
are considered to be accurate to within ±0.1°. 

The nose-boom installation for measuring the airspeed was calibrated 
by the NACA radar phototheodolite method. The Mach numbers presented are 
considered accurate to ±0.02 at speeds below about M = 0.90 and accurate 
to ±0.01 at speeds above M = 0.90. 

Instrument phase-lag corrections were applied to all data employed 
in the analysis. Also, position corrections were applied by time-vector 
methods of analysis to sideslip and to transverse acceleration data. 
Details of the application of the time-vector method axe considered in 
a later section of this paper. 


TESTS 


The test procedure for this investigation consisted of recording 
the airplane response to abrupt rudder pulses performed with other con- 
trols fixed. Attempts were made to maintain constant Mach number and 
altitude and to prevent inadvertent movement of the control surfaces 
during the transient portion of the maneuver. Such attempts were not 
always successful and required careful selection of usable portions of 
the flight record. Figures 6(a) and 6(b) present typical time histories. 
Small changes in altitude or Mach number did not appear to influence 
materially the results except in the region of the critical Mach number; 
however, moderate control movements in the transient portion of the 
maneuver influenced the analytical results. The most troublesome data 
resulted from maneuvers performed at high angles of attack or at other 
than 1 g . 

Maneuvers were performed at 1 g ±0.1g conditions for the four con- 
figurations at altitudes ranging from 38,000 to 4l,000 feet over a Mach 
number range of 0 . 73 to 1.35 • To extend the Mach number range of the 
tests to 1.48, maneuvers were performed following a pullout from a dive. 
These maneuvers were performed wjith configurations B, C, and D at 

35.000 ±3,000 feet over a load factor range of 1.2g to l-7g- 

To investigate the effects of angle of attack on the lateral sta- 
bility characteristics maneuvers were performed with configurations C 
and D during turns and pushovers at Mach numbers of 0.73 to l.l8 at 

40.000 ±2,000 feet and for configuration D at 30,000 ±2,000 feet. 
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Effects of angle of attack were also investigated over a Mach number 
range of 1.03 to 1.31 tor configurations A and B at altitudes from 39*500 
to 4l,500 feet with a load factor of 1.7g to 2.1g for configuration A, 
and an altitude range of 37*500 to 39*500 feet with a load factor of 2 . 1 g 
to 2.4g for configuration B. 


ACCURACY OF RESULTS 


In considering the probable errors in the analysis of the lateral 
characteristics of the airplane, attention must be given to instrument 
accuracy as well as to readability of the records, possible influences 
of variation in altitude and Mach number, influence of inadvertent move- 
ment of the controls, and accuracy of estimated derivatives. 

The readability of the records was a strong potential source of 
error. Since the ranges of the instruments and scale factors employed 
were governed by the roll-coupling investigations being conducted at the 
time, the deflections on the roll records were small in general and very 
small at Mach numbers in excess of about 1.3 • 

With all factors considered, the probable errors in the flight data 
employed in the determination of derivatives are estimated to be: 

Probable error, percent 


P: 

At M = O .75 0.5 

At M = 1.35 2 

T ± j 2 5 to 10 

o d 3 to 6 


M- 

i*r 


Subsonic region = 3-0 3 

I'M 

Supersonic region -pM- = 1.6 6 

m 


The probable errors in the lateral stability derivatives obtained 
from flight data are dependent on the degree of error in the estimated 
values of Cj. and in the moments of inertia, and in the direction 
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of the principal axis, as well as the errors presented in the previous 
paragraph. The probable errors in the derivatives, exclusive of the 
influence of errors in the estimation of C n and C z , which will be 

.fcr 

discussed later, are believed to be: 

Derivatives Probable error, percent 


' n P 


Cl 


3’ 

From M = O.7O to 1.30 
At M = 1.47 


'h 


From M = O.7O to 1.30 
At M = l.li-7 


'n T 


- C 


n fi) 


3 to 5 
10 


10 

50 


10 

20 

15 to 20 


PROCEDURES OF ANALYSIS 


The original wing area and wing span were employed for all config- 
urations in analyzing the flight data for the lateral stability deriv- 
atives . To convert the derivatives of the extended wing configuration 
(configuration D) to the actual wing area and wing-span bases, the 
Cy R derivative should be multiplied by 0.975; and C i a by 0.925; 


3 

and 


C 7 and C n by O.876. 


l 3 


Inasmuch as the maximum sideslip and roll angles of the tests were 
of the order of 2° and 10°, respectively, and since no significant non- 
linear or cross-coupling influences were noted, the following linearized, 
small disturbance forms of the lateral equations of motion of the air- 
plane appeared applicable to the analysis of the data: 

Wa-t = m(v + rV - gcp) (l) 


Wat ~ ( C YpP + C Y p P' + Cy r r ' + Cy^P'^qS (2) 

%P " X XZ* = (%3 + C Zp p' + C Zr r’ + C z .fi')qSb ( 3 ) 

1 7? r ~ X XZP = ( C n p P + c n p P' + c n r r ' + c n^'j<lSb (4) 
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Three methods were originally considered for the determination of 
the lateral stability derivatives. The frequency -response method of 
reference 8 was highly desirable because of the number of derivatives 
which can be determined from it; however, because of the time factor 
and some doubt as to the validity of the results which would be obtained 
using the available flight data, it was decided not to employ this method. 
The method of reference 9 is a time-vector approach to the solution of 
the derivatives; however, it is a tabular procedure employing successive 
approximations and therefore is not as desirable as the relatively rapid 
graphical time-vector method of analysis explained in references 10 to 12. 


for 


The graphical time-vector method of references 10 to 12 was employed 
the determination of Cy^, C^, ^C nr - C^ , C z ^, and C Zp . The 


required precision of phase-angle data precluded the possibility of 
reliable values of ^0y r - Cy^j or Cy^; therefore it was decided, on a 

selective basis, to employ estimated values of Cy p and to ignore 


( 


- Cy^ in the solution. The values of C ^ and C n ^ which were 

required for the time-vector solution of the other derivatives were 
obtained from theoretical estimates. 


Application of the Time-Vector Method of Analysis 

No attempt is made in this paper to present the detailed mathematical 
aspects of the fundamental time-vector properties inasmuch as reference 10 
accomplishes this quite thoroughly. Suffice it to say that the time 
invariance of the phase relationships and amplitudes relative to each 
other permits the representation of any one of the linearized equations 
of motion by vectors. In the four lateral-directional equations three 
degrees of freedom are involved in each equation; namely, sideslip, roll, 
and yaw, each with the same frequency and damping characteristics. The 
amplitudes of the various degrees of freedom in each of the lateral- 
directional equations have the same shrinkage rate and the phase angles 
remain constant; thus for vector representation, the various amplitudes 
and phase relations are time invariant. 

The vector properties described in the preceding paragraph, plus the 
requirement that the vector polygon representing any one equation must 
close, makes possible the determination of two unknowns in any one equa- 
tion. Inasmuch as it is desired to determine the stability derivatives 
from flight data, it will be convenient to introduce new notations for the 
stability equations and to establish the equations in the form of ampli- 
tude ratios. All equations in this paper having absolute value notations 
will be considered to represent vector equations. Hence 
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(6) 


(5) 




(8) 


The derivatives with respect to r and (3 have been combined in 
equations ( 6 ), ( 7 ), and ( 8 ). This was done because^ |r| is similar to 
|( 3 | and is approximately l 80 ° out of phase with |p| . 

The amplitude ratio representation is convenient, inasmuch as it 
simplifies flight-data reduction and enables a more direct determination 
of some of the derivatives. 

The period of oscillation P is determined directly from the tran- 
sient portion of the flight record. To determine the indicated phase 
angles, the measured time differences of the different peaks of the 
various degrees of freedom were averaged and the simple expression 


was utilized. 

To determine indicated amplitude ratios relative to the body axes, 
the envelopes of the transient oscillation records are plotted on a semi- 
logarithmic plot, such as figure 7 , as a function of time. The linearity 
of the curves indicates that the linearized differential equation is 
applicable. The indicated amplitude ratios are then calculated as shown 
in figure 7 - 

The values of T]_/2 are determined as indicated in figure 7* The 
damping angle % is related to the damping ratio £ by one of the 
relations 


* = ± ^( 360 °) 


(9) 
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= tan - ^- £ 

l/rrp 


= sin" 1 ^ 


= COS 


_1 l/l - ( c 


(10) 


The determination of or £ requires the use of the logarithmic 

decrement which may he expressed in various ways such as 


S = ^ _ _ = 2rt tan 


h-t‘ 


and 


5 =2^23? 


T l/2 

Equating (ll) and (12) and transposing 


( 11 ) 


( 12 ) 


tan % = °; - 69 - 3 — = 0.1102 P 
2rt T l/2 


T l/2 


Since cd^ = o^^l - , then from (10) 


(13) 


con = — — 


2rt 


COS <J>£ p cos % 


(14) 


Correction of Indicated Amplitude Ratios and Phase Angles 

Amplitude ratios are subject to corrections for dynamic magnifica- 
tion, instrument location, and reorientation when the data are to be 
considered relative to axes other than the body axes about which the 
instruments are oriented. Phase angles are subject to corrections for 
phase lag in the electronic system, phase error caused by instrument 
location away from the center of gravity of the airplane, and reorien- 
tation when the data are to be considered relative to axes other than 
the body axes. The phase lag of the instruments was based on the 
relation 
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Phase lag = -tan 


-1 


2 W% 



deg 


(15) 


where 


a) frequency of the airplane, radians/sec 

o) ni undamped natural frequency of instrument, radians/sec 

The indicated phase angles were corrected for the difference in the phase 
lag of the instruments involved. 

The transverse accelerometer and (3-vane records were subject to phase 
error caused by instrument location. These phase errors were determined 
in a vector approach to the correction of indicated amplitude ratios to 
true amplitude ratios. 

The amplitude ratios will, in the case of carefully selected instru- 
ments, be subject to negligible dynamic amplification error. In the case 
of properly oriented gyroscopic instruments, the position error will be 
negligible and location error is not a factor to be considered. In the 
case of the transverse accelerometer and the (3 -vane instruments, location 
error may be important as on the present airplane. 

To correct the (3 record for vane location, the following expression 
was employed 


3 = Pi 


x p r b , z pPb 

V V 


In terras of vector notations 

iiL - , iJM 
N ' " v |»il v N 


( 16 ) 


( 17 ) 


The graphical time-vector solution of ( 17 ) is shown in figure 8(a). 
The solution is obtained by first drawing in the direction of the vectors 
relative to (3 using the indicated phase angles corrected for instrument 
phase lag. The numerical values of the terms in ( 17 ) are then drawn in 

M 


as vector quantities . The magnitude of the 


IN 


vector and its direction 


represent the magnification factor by which all the amplitude ratios taken 
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with respect to 0 should be divided to correct for 0 position error, 
and the phase angle error of all of the phase angles taken with respect 
to 0 . 


To correct the indicated transverse acceleration for position error 
the following expression was used 


a t = a t n + 


z a^b x a^_ r b 


I g 


g 


In terms of vector notations 


1 = foil , Za t N Xa t It 


where 


IPI |Pl § IPI g |P| 
I Pbl _ Kl 


lei | el 


“n 


and 


IPI I P I 


cpn 


(18) 


(19) 


The time-vector solution of ( 19 ) is shown in figure 8(b). The solution 
is approached by first drawing in the directions of the p^ and r^ vec- 
tors . The directions of the acceleration vectors p^ and are then 

located (90° + ahead of the velocity vectors. The remainder of the 

solution is as shown in figure 8(a). It can be seen that the position 
correction is significant in the illustration representing actual test 
data of an average condition. 


Determination of the Lateral Stability Derivatives 

The stability derivatives being determined and the equations of 
motion employed in the present analysis are referenced to the stability 
axes. Inasmuch as the flight data are referenced to the body axes, it 
is necessary to transfer the flight data from the body axes to the sta- 
bility axes. Had the stability derivatives sought been referenced to the 
body axes, then the stability equations referenced to the body axes would 
have been employed. In the appendix are presented the equations for the 

transfer of the amplitude ratio j^-j and the phase angle $p r from the 
body axes to the stability axes. The appendix also outlines the vectorial 


CONFIDENTIAL 


18 


CONFIDENTIAL 


NACA RM H56C20 


procedure employed in obtaining refined magnitudes of 


Hi 

IP 


and 


as 


well as phase angles with the aid of the transverse acceleration equa- 
tion (5)- Figure 8(c) illustrates the final vector solution for these 
magnitudes . 

After correcting the data for various sources of error and trans- 
ferring the data to the stability axes, as shown in the appendix, it is 
a simple matter to proceed with the determination of the derivatives . 

Since the positions of the displacement, velocity, and acceleration vec- 
tors are known, the three lateral-stability equations ((6) to (8)) may 
be applied to the solution of the lateral derivatives. 

Figure 8(d) shows the vectorial solution for Cy^ and ^Cy^ - Cy^ * 
Although a solution for ^Cy r - Cy^ is shown in the figure, this deriv- 
ative was not included in the results of the analysis because of the lack 

of the required precision of the value of $ a ^ which would be needed 

t 

to obtain a fairly reliable first approximation of this derivative. 

Figure 8(e) shows the vector diagram for the determination of Cn^ 

and ^C nr - C n ^ . No attempt was made to determine C n .p in place of one 

of the other two derivatives, since some preliminary work appeared to 
indicate there would be no advantage in doing this. The section entitled 
"Discussion" in this paper considers sensitivity of some of the derivatives 
to experimental errors as well as limitations in the application of the 
time-vector method of analysis. Figure 8(f) shows the vector solution 
for C jp and C^. In this solution was neglected since 


Cl r “ (< 


- %)■ 


Estimated Derivatives 


For the static derivatives 


C Y P 


and 


tp, “ 1U - % 

estimates based on wind-tunnel data were obtained from manufacturer’s 
estimates. Wing contributions to the dynamic stability derivatives were 
estimated from the methods of references 13 to l8. The vertical-tail 
contributions to the static and dynamic stability derivatives were calcu- 
lated, using the method of reference 19 and calculated lift-curve slopes 
(refs. 14, 17, and 20 ). 


available tail-off 


Manufacturer's estimates were utilized (fig. 9) for the effect of 
torsional flexibility on vertical-tail contributions to C v . . Cn., C 


and Cy for tail C. 


' y 3 * 

The flexibility corrections were also applied to 


n r > 
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the calculation of 


C Z R - 


Figure 10 shows the estimated change in the 


derivatives caused hy vertical-tail torsional flexibility in configura- 
tions C and D. 


Two sidewash influences were considered to be acting on the vertical 
tail. The sidewash factor caused by roll Op was, on the basis of ref- 
erence 19, estimated to be 0.25- The sidewash factor caused by side- 
slip Op has been shown in references 21 and 22 to be a function of wing 

position and influences the values of Cy^ > c n p > (Cn r “ c np) • 

From reference 22, Op was estimated to vary somewhat linearly for the 

angle -of -attack range of flight tests and was considered to vary from 
approximately O.O 7 at a = 0° to 0.11 at a = 6 °. 


The side force at the engine duct inlet resulting from the momentum 
change caused hy bending of the intake air to flow along the duct axis 
was also taken into account in calculating Cy^ and of configur- 

ations C and D. No air-intake effects are considered when calculated 
derivatives on a rigid tail basis are presented. The air-intake effects 
were approximated by the equations 


(^)a ' T 

(20) 


(21) 


Figure 
engine 
but ion 


11 shows the estimated weight rate of air required by the jet 
to maintain cruising speed. Figure 12 shows the estimated contri- 


of the intake air of the jet engine to Cy^ and C 




PRESENTATION OF RESULTS 


Following is a summary of the figures presenting the results of 
this investigation: 


Limitations of the Time-Vector Method 



Figure 

Influence of and C ^ 

13 

Influence of ±l-percent change in Cnp 

ik 

Influence of ±5-percent change in $p r 

15 

Influence of ± 0 . 5 ° change in 0^ 

16 
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Static and Dynamic Characteristics 



• 


Figure 


Flight data 
results 

Configuration 

Trim, 

a 

Period 

and 

damping 

Amplitude 
ratios and 
phase angles 

Static and 
dynamic lateral 
derivatives 

A 

17 

18 

19 

20 

B 

17 

21 

22 

23 

C 

17 

24 

25 

26 

D 

Comparison of A, B, 

17 

27 

28 

29 

C, and D 

Comparison of C and D 
with theory and 


30 

31 

32 

wind tunnel 

Angle-of-attack effects 
at M = O.83 and 




33 

M = 1.12 

— 

34 

34 

35 


The data for configuration A, shown in figures l8 to 20 , are meager 
in the subsonic region and most are subject to inadvertent control move- 
ments which, although not affecting the periods (fig. 18(a)) appreciably, 
do affect the damping (fig. 18(b)) and the phase angles (fig. 19(b)) so 
that no attempt was made to analyze these data for the 40 , 000 -foot con- 
dition. The three test points at M = O.7I constitute the only reliable 
damping characteristic points in the subsonic region and, as a result, 
the amplitude ratio curves of figure 19 indicate approximate values only. 
Despite the lack of sufficient subsonic data, the experimental stability 
derivative characteristics shown in figure 20 are considered to be reli- 
able within the accuracy indicated previously. 

Although period and damping curves are shown in figure l 8 for a 
load factor of 1.8 at 40,000 feet, the amplitude ratios and phase angles 
for this condition were not sufficiently well defined to obtain derivatives. 

The results of the analysis for configurations B, C, and D (figs. 21 
to 29) are based on the availability of a larger amount of pulse data for 
each configuration. The data for configuration C were sufficient to 
define characteristic curves for trim level flight at 31,000 feet from 
M = 0.77 to M = 1.0 as well as for trim level flight at 40,000 feet 
(figs. 2 k to 26). 
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DISCUSSION 

Limitations in the Application of the Time-Vector Method of Analysis 


Although the time-vector method of analysis requires the simplest 
of equipment in its application and is capable of providing good results, 
it does have definite limitations. In considering the limitations, it 
is presumed that flight records have clear-, sharply defined traces, and 
that the ordinates have ample magnification in relation to the period 
scale to produce well-defined peaks in the oscillations. 

One of the limitations in the application of the method is the 
inability to work with records of heavily damped airplanes without 
resorting to other methods of analysis, such as frequency -response anal- 
ysis, to obtain amplitude ratios, phase angles, and angular frequency 
of the motion; and the use of template aids or analogs to determine the 
damping ratio. When the damping ratio £ exceeds 0.2, the accuracy of 
defining Tq /2 ^> e Si ns to decrease. When £ exceeds approximately O.jJO, 

it is somewhat difficult to determine the period accurately and the ^ 1/2 

values become increasingly doubtful. Also when £ exceeds 0.4, relia- 
bility of P and T -^2 becomes poor. 

For controls-fixed conditions, the method depends on the analysis 
of the transient portion of an oscillatory motion. Any inadvertent 
application of a forcing function during this transient oscillatory 
motion, even though it may be small, will tend to influence the results. 
In instances where the forcing function is deliberate and is of a pure 
sinusoidal nature, the time-vector method is applicable providing the 
Cyg , , and C^g derivatives are available. 


A third limitation of the time-vector method lies in the fact that 
only two of the three derivatives in each of the lateral equations may 
be determined by means of the vector diagram. 


and 


In the case of transverse equation (6), the secondary terms 




are generally neglected and the result is 



This simplified expression for Cy 


JL ° 01 

provides answers which are high} 


(22) 


however, the error probably does not exceed 4 percent. The principal 
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difficulties in obtaining refined values of Cy^ have been in the read- 
ability of the records and the phase lag error of the vane itself. It 
should be recognized that the unknown phase lag of the vane would enter 
into the problem and affect the answers for Cy^, regardless of the 

method of analysis employed. 

In the case of the rolling-moment equation (eq. (8)) it was deemed 
advisable to estimate the values of Cj r and to obtain and 


from the vector diagram. The Cj 


r IPI 


vector is relatively small, 


especially at high Mach numbers, and a normal error of ±5° in would 

result in no accuracy in attempting to determine Cj . 


A limited investigation was made of the sensitivity of the deter- 
mined derivatives to variations of the assumed values of C n p and C i r , 

to ±5° errors in $p r > to ±0.5° error in damping angle, and to a ±1 per- 
cent change in C n p. This investigation was considered for configura- 
tion C at M = 0.80 and M = 1.20 at 4-0,000 feet. 

As shown in figure 13, has a relatively small effect on C n ^ 

and a fairly large effect on ^C nr - C n pj • The effect on ^C nr - C n p^ was 

of the order of 0.029 and 0.017 units per 0.01 unit change in C n at 

M = 0.80 and M = 1.20, respectively. The results of figure 13 show 
that is affected less than 0.0004 units per 0.01 unit change in 

C and that C ^ is affected 0.0037 and 0.0052 units at M = 0.80 

and 1.20, respectively, per 0.01 unit change in C j . 


Normally, in dealing with the yawing-moment equation (eq. (7)), 
attempts are made to determine the ^C nr - derivatives from vector 

diagrams. Thus, either C n p or C n p must be obtained by other means 

to permit completion of the solution. In the present paper a theoret- 
ical estimate of C n p was made and used to obtain both ^C nr - C n ^ and 

. Inasmuch as there is usually some question of the accuracy of C n ^ 
estimations, some vector solutions of C n p and ^C nr - C n ^ were obtained 
using as determined from 


"np 



_ 1 _ 

qSb 

P 2 


JL + 1 f 0.693 ^ 


4it 2 \ T l/2 


- Cl 


P I 


^-xz 

X 


(23) 
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In considering the possibility of employing calculated values of C 


to determine C 




and 


“P 


(^n r “ ^nj ) y values of C n were obtained from 


equation (23) using faired flight data and vector solutions of C 7 . 

The influence of a ±l-percent change in on the lateral stability 

derivatives C n p and (C nr - C n ^ is shown in figure 14. On the basis 
of the vector diagrams shown in figure 14, the influence would be appre- 


ciable. In view of the influence of small errors in 
effects of the possible magnitudes of errors in $ 


' n p 


plus the 


p r and Od it was 
decided that, insofar as the present analysis is concerned, it would be 
better to employ theoretically estimated values of C 


n T 


Although the errors in the phase angle are believed to be 

generally within +3° in the present paper, the error may approach ±5°. 
Figure 15 shows that a ±5° error had negligible influence on C n p, had 

a moderate effect on Cjp, pronounced effects on C ^ and - C n ^) 

A study of the vector diagrams in figure 15 will show that a decrease in 
the magnitude of will reduce the influence of phase angle error 

in C n p and ^Cn r - Cn^), but will also increase the influence on Cj 
and C i . When the phase lag of the p vector decreases toward 

90 °, the influence of the t5° error on both C i and C i increases. 


A ±0-5° error in the damping angle 0^ showed small to moderate 
influence on ^C nr - C n ^, as shown in figure 16. In instances where 
^Cn r - Cnp) would be of the order of -0.10, the error would be pronounced. 


Comparison of the Four Configurations 

A summary of the results of the analysis of the flight data of the 
four configurations to show the influences of the various modifications 
on the stability characteristics is presented in figures 30 to 32. 

The period characteristics (fig. 30(a)) show an appreciable decrease 
in the period when the original tail of configuration A was extended to 
form configuration B. Replacement of the extended tail by the large tail 
to form configuration C showed a moderate decrease in the period over the 
Mach number range. The extension of the wing in configuration C to form 
configuration D had a small uncertain effect on the period. The effects 
of the various modifications on the period characteristics are reflected 
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in the characteristics of the directional stability derivative Cnp as 
seen in figure 32. 


Although the damping characteristics (fig. 30(b)) show that 

was decreased in the subsonic range and increased in most of the super- 
sonic range with each increase in tail size, the damping ratio £ showed 
a decrease with each increase in tail size throughout the entire Mach 
number range investigated. The addition of wing-tip extensions, to form 
configuration D, appears to have negligible effect on both T £ 

from M = O.7I to about M = 0-9- Between M = 0.9 and M = 1.37; the 
addition of the wing tips appears to increase the damping. 


The influence of the increase in vertical-tail sizes and the addi- 
tion of the wing-tip extensions on the damping parameters T-^y^ and £ 

is perhaps most effectively expressed in terms of derivatives as shown 
by approximate relationships based on the analytical expressions of 
reference 23- Although not exact, the following relationships, applicable 
to low angle-of-attack conditions, appear generally adequate for quali- 
tative purposes: 



1.386 


gS 

V 


2I7 (° n r 




/qSI z 1 


C »r - S) + * °*fl) 

V 4b 



(24) 


(25) 


A study of equations (24) and (25) indicates that at any one value 
of dynamic pressure, T^y 2 is 

ratio £ is dominated by both 

derivative characteristics (fig. 32(b)) show qualitative trends with 
configuration which are, in general, compatible with the T^y 2 trends 

shown in figure 30 . 


dominated by and the damping 

^Cn r - ^ n p) * •^ ie (^n r ~ 0 n ^ 


The pertinent amplitude ratios are shown in figure 31- The character- 

I a + I 

istics curves of L , - V - for the various configurations are somewhat erratic 

im 

relative to each other because of the poor readability of the transverse 
acceleration flight records and the p-vane errors discussed previously. 
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In the subsonic range the large values of 


of configuration A 


were markedly reduced by the various tail and wing modifications. In 
the supersonic region configuration A showed the lowest j^-j magnitudes 
and an increase in vertical-tail area increased the ratio; however, the 
addition of extended tips decreased |^-j- slightly. 


The phase angle 0^^ did not appear to be influenced in the sub- 
sonic range by the range of vertical-tail sizes covered, but the addi- 
tion of the wing extensions had a more significant influence on the phase 
angle (fig. 3l). In the supersonic region configuration A showed less 
lag in phase angle than did configurations B and C, which had practically 
identical phase-angle characteristics. Extension of the wing tips tended 
to decrease the lag. 


In the Mach number range beyond M = 1.2 or 1.25 the amplitude 
ratio and phase-angle characteristics appear, in general, to be changing 
at an increasing rate. These changes in characteristic trend are reflected 
in the derivative characteristics shown in figure 32. 

Figure 32 shows that an increase in both vertical-tail size and 
aspect ratio had desirable influences in the trim level-flight static 
derivative characteristics. Configuration C had practically double the 
directional stability of configuration A at M = O.7, and approximately 
a 70-percent increase throughout the supersonic range. The influence 
of the different vertical tails on the directional stability has been 
reported previously in reference 2 relative to the body axis. When the 
C R p curve of this paper for configuration C and 40, 000-foot altitude 

was transferred to body axes and compared with reference 2, excellent 
agreement was evident over the entire Mach number range. 


The effective dihedral 


C L: 


was also subject to substantial increases 


with each increase in vertical-tail size. Wing-tip extensions had negli- 
gible effect. The rather sharp reduction in the negative value of 

in the vicinity of the critical Mach number of about O.96 is caused by 
the tail-off characteristics of the airplane. The deterioration of 
effective dihedral with increasing Mach number from M « 1.23, when con- 
sidered in conjunction with the break in the curves for configura- 

P 

tions C and D, tends to indicate the possibility of shock wave and flow 
interference near the tips of the wings which influences the lift distri- 
bution across the span of the wing. Such an influence would reduce the 
effective dihedral C ^ which tends to become negative in the region of 

M = 1.38 to M = 1.47, depending on the configuration. 
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The major influence of wing -tip extensions (configuration D) appeared 
to be on the Ci derivative. In the supersonic range between M = 1.05 

and 1.30, there appears to be a fairly large increase in the negative 
value of the damping-in-roll derivative Cj . This indicated increase 

is based on the original wing area and span. When based on the actual 
wing area and span, the damping in roll for configuration D is larger 
than for configuration C up to M « 1 . 31 - 

In view of the difficulty in obtaining ^C nr - ^ n p) derivatives arLd 

in an effort to check roughly the magnitudes of the - ^ n p) deriv- 

atives as obtained by the time-vector method, equation ( 24 ) was trans- 
posed to the following form to solve for (C nr - C n -) : 



c \ ~ 2I Z /1.386V 
V " b 2 [qST l/2 



(26) 


Utilizing the Tp/2 311(1 C Yp values for configuration D (figs. 30 
and 32), (C nr - C n .) was computed by using equation (26) and compared 
as shown in the following tabulation, with the values of ^C nr - Cn^ as 
determined by the time-vector method. 


Mach number 

0.90 

1.05 

1.15 

1.25 

(Cn r - Cn^) (eq. (26)) 

- 0.34 

- 0.257 

- 0.239 

-0.262 

(Cn r - by vectors 

- 0.29 

-0.200 

-O.205 

-0.210 


Inasmuch as equation (2 6) is approximate and tends to provide 
(Cn r ~ c np) values which are high (especially at higher angles of attack), 

it appears that the vector solutions for ^C nr - Cn p) are vi thi n reasonable 

limits and a rough insight as to the influence of the various configura- 
tions may be justified. 


The negative magnitude of (Cn^ 



appears to increase with 


increase in vertical -tail size in the subsonic range. Supersonically 
there appears to be a decrease in negative magnitude with increase in 
vertical-tail size. The addition of wing-tip extensions decreased the 
negative magnitude of ^C nr - 1° some extent; supersonically the 

influence appears to be negligible. 
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Comparison With Calculated Characteristics and 
Wind-Tunnel Data 

Two sets of calculated characteristics curves are shown in figure 33* 
The results show that air-intake effects and torsional flexibility of 
the tail have a pronounced influence on the calculated stability 
characteristics . 


Beyond M = 1.25 all the flight-determined derivatives except 
( c n r - c np) experience a deteriorating break in magnitude characteristics. 

The calculated Cy 


n o characteristics indicate this break clearly; 


calculated C 


U 


and C 
3 “3 

characteristics show only slight but similar trends 


starting at M = 1.15, calculated characteristics indicate that 

damping in roll begins to deteriorate in the vicinity of M = 1.35* 

Inasmuch as is practically dependent on wing alone, the break 

in the C ^ curve not accounted for by calculated values of this deriv- 
ative appears to indicate, as mentioned in the previous section, the 
possibility of some shock wave and flow interference near the tip of the 
wings of both configurations which influences the lift distribution 
across the span of the wing. Such an influence would reduce the effec- 
tive dihedral which tends to become negative at a Mach number of 

approximately 1.47* 


A comparison of the calculated derivatives with flight results 
showed fair to good agreement in the subsonic region for all derivatives 
except ^C nr - C n ^ . The calculated values of ^C nr - C n ^ , similar to 

the low-speed wind-tunnel values, were much lower than flight results. 

Unpublished wind-tunnel static-stability data for M = 1.4-1 were 
corrected for vertical-tail flexibility and air-intake effects of the 
jet engine and are plotted in figure 33. These modified wind-tunnel data 
show good agreement with the flight-determined trend of Cn^ and C j . 

It is difficult to compare the low-speed wind-tunnel data with the 
subsonic flight results (fig. 33) because of the large Mach number differ- 
ence. As will be pointed out in the following section, the variation 
of Cy^, C n ^, and Cj with angle of attack shown by wind-tunnel data 

is the opposite of trends shown by flight results; however, it appears 
that the magnitudes of and C z from flight and wind-tunnel data 

tend to agree. 


CONFIDENTIAL 


28 


CONFIDENTIAL 


NACA RM H56C20 


Variation of Stability Characteristics With Angle of Attack 

Although available flight data permitted the presentation of con- 
stant load, factor curves for several of the lateral characteristics in 
the supersonic region for configurations A and B (figs. l8 to 21 ), and 
data were also available for the presentation of the lateral character- 
istics in the subsonic region for altitude effects for configuration C 
(figs. 24 to 26 ), no attempt is made in this paper to discuss the results 
inasmuch as other flight data provided a more detailed insight into the 
variation of the lateral stability characteristics with angle of attack 
at Mach numbers of 0.8l and 1.14 . 


The variation of the lateral stability characteristics with angle 
of attack at M = O .83 at altitudes of 40,000 and 31,000 feet for con- 
figurations C and D, and at M = 1.14 at an altitude of 40,000 feet for 
configuration C are shown in figures 34 and 35- Also shown in figure 35 
are the variations of C n p, C^, (^ n r “ ^ n p) a^gde of 

attack as obtained from reference 1 for a Mach number of 0.13. 


As shown in figure 34, flight data indicate a decrease in period 
with increasing angle of attack regardless of the Mach number or altitude. 
The damping characteristics improve with both increasing angle of attack 
and decreasing altitude. 

The amplitude ratio j^j- and the phase lag of fc^p increase with 

angle of attack. Increasing angle of attack tends to place the roll 
and sideslip displacements in phase. This tendency, plus the increase 
in roll angle per unit sideslip angle, tends to accentuate Dutch roll 
tendencies of the airplane. 


Figure 35 shows the trends of the variation of the derivatives with 
angle of attack. The Cy derivative is not included because the scatter 


of the flight results precluded the possibility of presenting a definite 
trend of Cy^ variation with angle of attack at constant Mach number. 


Although low-speed wind-tunnel data from reference 1 are shown for com- 
parison with flight results at M = O .83 and M = 1.14, a direct compar- 
ison for the same Mach number conditions is difficult because of the 
presence of automatic leading-edge slats on the airplane and the large 
difference in Mach number which would make extrapolation unreliable. 


Flight results indicate an increase in directional stability and 
effective dihedral with increasing angle of attack. 


The damping-in-roll derivative 
magnitude at an angle of attack of about 


appears to attain its maximum 
3°. 
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CONCLUSIONS 


From the analysis of flight data obtained for four configurations 
of a swept-wing fighter-type airplane over the Mach number range from 

O .7 to 1.14-8 the following conclusions have been reached: 

1. The time-vector method of analysis is capable of producing good 
values of the lateral derivatives Cy^, C n p, C^, anc ^ ^Zp providing 

the damping ratio is less than approximately 0.3- Reliable values of 
lateral derivatives ^C n ^ - C n ^ are difficult to determine because of 

the sensitivity of this quantity to other factors. 

2. The expected effects of increasing vertical-tail size, resulting 

in increased magnitudes of C^, C^, and C^, were realized. 

3. The addition of wing-tip extensions had small effects, except 
for a fairly large increase in the magnitude of the damping-in-roll 
derivative C 7 . 

4. Theoretically calculated derivatives showed fair to good agree- 

ment with flight results in the subsonic range with the exception of 
high angle-of -attack values of derivatives. Wind-tunnel data 

for the static derivatives for a Mach number of 1.4l, when corrected for 
torsional flexibility and air-intake effects of the jet engine, showed 
good agreement with flight results. 


C 


v 


5« The experimental rate of decrease in the magnitudes of Crip* 
and Cjp with Mach number at Mach numbers greater than 1.25 vas 


larger than estimated. This increased rate of decrease in magnitudes 
appears to be the result of possible shock wave and flow interference 
at the wing tips . 


High-Speed Flight Station, 

National Advisory Committee for Aeronautics, 
Edwards, Calif., March 9 , 1956. 
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APPENDIX 

TRANSFER OF, AMPLITUDE RATIOS AND PHASE ANGLES 
FROM BODY AXES TO STABILITY AXES 


The transfer of the amplitude ratios and phase angles from body axes 
to stability axes was accomplished by the use of equations for the trans- 
fer of !M and <D to the stability axes, and the use of the vector 

I ^b b 


method to obtain j-^-, $^0, and 


To transfer 


lei 

|p- 




and 3>p^ r ^ to the stability axes, the following 


equations from reference 9 were employed 


M = 

Irl 



(Al) 


and 


$ = tan 

pr 


-1 


^b 


sin <J>- 


p b r b 


cos 0. 


p b r b 


+ tan a 


tan' 


■M 


p b 

— - sin <tv, „ tan a 
r ul Pb r b 


1 - H cos .. 


p b r b 


tan a 


The amplitude ratios 


and 


M 

ip 


(A2) 


and the phase angles $ r p and 


$^p were obtained vectori ally with the aid of the transverse acceleration 

equation (5) as shown in figure 8(c). In approaching this analysis the 
directions of the £4. and p vectors are drawn as shown and, since nei- 


ther 


m 


nor the direction of the cp vector is known, first approxi- 


mations are made for these quantities as follows: 
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and 


where 


M ^ |p I 

fPf ~ 

V ~ ^cp\|r + 180 ° 



180 ° 


M = M 

|r| ji|rj 


and p 


$ 


(A3) 

(A4) 

(A5) 


As a result of these first approximations , the vector diagram is 
drawn as shown. The closing vector 2t j~- determines the first approxi 

mation of the direction of the r vector from which the second approxi- 
mation of the \| t vector direction is determined to obtain the second 
approximation of the cp vector direction using 4 pr . 


To obtain the second approximation of divide 2r -~j- by 

2t to obtain the first approximation of -^j- and multiply this 


value of 


* 


Using the second approximation of 


tor, the second approximation of C 


io 


q> 


<p 


• and direction of the \|r vec 
is determined and redrawn on 


the vector diagram to obtain a new value of 2t j-£-j- and direction of the 
r vector. This second approximation of 2t is now used to obtain 


a second approximation of 


* 


I cp I 

and a third approximation of 


It has been found that the direction of the r vector and the 

magnitudes of |^-j- and determined by carrying the successive 

approximations thus far are quite close to the values which would be 
obtained had the successive approximation procedure continued to complete 
convergence. 


Having determined 


♦ I 


1 9 


W IP 


, and the direction of the r vector, 


it is a simple matter to obtain the phase angles 4^ and 4^. 
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TABLE I.- PHYSICAL CHARACTERISTICS OF AIRPLANE 


Basic wing 


Wing: 

Airfoil section 

Total area (including aileron and 83.84 sq ft covered By fuselage), sq’ft ! ! 

Span, ft ' 

Mean aerodynamic chord, ft 

Root chord, ft * 

Tip chord, ft 

Taper ratio 

Aspect ratio 

Sweep at 0.25 chord line, deg 

Incidence, deg 

Dihedral, deg 

Geometric twist, deg 

Aileron: 

Area rearward of hinge line (each), sq ft 

Span at hinge line (each), ft !.!!!!!!! 

Chord rearward of hinge line, percent wing chord 

Travel (each), deg 

Leading-edge slat: 

Span, equivalent, ft 

Segments 

Spanwlse location, inboard end, percent wing semi span 

Spanwlse location, outboard end, percent wing semispan 

Ratio of slat chord to wing chord (parallel to fuselage reference 

line), percent 

Rotation, maximum, deg 

Horizontal tail: 

Airfoil section 

Total area (including 31 . 65 sq ft covered by fuselage), sq ft 

Span, ft 

Mean aerodynamic chord, ft 

Root chord, ft 

Tip chord, ft 

Taper ratio 

Aspect ratio 

Sweep at 0.25 chord line, deg 

Dihedral, deg [ ’ ’ ’ 

Travel, leading edge up, deg .'!!!!!.*!! 

Travel, leading edge down, deg !!!!!!.* 

Irreversible hydraulic boost and artificial feel 


NACA 64A007 
376.02 

36.58 

11.33 

15.86 

4.76 

0.30 

3.56 

^3 

0 

0 

0 


19.32 

7.81 

25 

±15 

12.71 

5 

24.6 

94.1 

20 

15 


Basic wing plus 
wing-tip extensions 


NACA 64A007 
385.21 
38.58 
11.16 
15.86 
4.15 
0.262 
3.86 

^3 

0 

0 

0 

19.32 

7.81 

25 

*15 

12.71 

3 

23.3 

89.2 

20 

15 


NACA 65A003.5 
98.86 
18.72 
5.83 
8.14 
2.46 
0.30 
3 . 5 ^ 
^5 
0 
3 

25 
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TABLE I.- PHYSICAL CHARACTERISTICS OF AIRPLANE - Concluded 


A 


B 


C 


On 


Vertical tail: 

Airfoil section 

Area (excluding dorsal fin and area blanketed by fuselage), sq ft . 
Area blanketed by fuselage (area between fuselage contour line and 
line parallel to fuselage reference line through intersections of 
leading edge of vertical tail and fuselage contour line) .... 

Span (unblanketed), ft 

Mean aerodynamic chord, ft 

Root chord, ft 

Tip chord, ft : 

Taper ratio 

Aspect ratio 

Sweep at 0.25 chord line, deg 

Rudder: 

Area, rearward of hinge line, sq ft 

Span at hinge line, ft 

Root chord, ft 

Tip chord, ft 

Travel, deg 

Spanwise location, inboard end, percent vertical tail span . . . 
Spanwise location, outboard end, percent vertical tail span . . . 

Chord, percent vertical tail chord 

Aerodynamic balance 


NACA 65A005.5 

NACA 65A005.5 

NACA 65A005.5 

33.5 

37-3 

42.7 

2.11 

2.11 

2.45 

6.14 

7-^5 

7-93 

5.85 

5.51 

5.90 

7-75 

7*75 

8.28 

3.32 

2.52 

2.49 

0.428 

0.501 

0.501 

1.13 

1.49 

1.49 

45 

45 

45 

6.5 

6.5 

6.5 

3.33 

3.33 

3-33 

2.27 

2.27 

2.27 

1.50 

1.50 

1.50 

t20 

±20 

±20 

4.5 

3.7 

3.1 

58.2 

48.0 

44.8 

50.0 

50.0 

28.4 

Overhanging, 

Overhanging, 

Overhanging, 

unsealed 

unsealed 

unsealed 


Fuselage: 

Length (afterburner nozzle closed), ft . 

Maximum width, ft 

Maximum depth over canopy, ft 

Side area (total), sq ft 

Fineness ratio (afterburner nozzle closed) 


45-64 

5.58 

6.57 

250.92 

7.86 


Speed brake: 

Surface area, sq ft . . 
Maximum deflection, deg 


14.14 

50 


Power plant: 

Turbojet engine 

Thrust (guarantee sea level), afterburner, lb 

Military, lb 

Normal, lb 

Airplane weight, lb: 

Basic (without fuel, oil, water, pilot) 

Total (full fuel, oil, water, pilot) 

Center-of -gravity luoation, percent c: 

Total weight - gear down 

Total weight - gear up 

For moments of inertia and inclination of principal axis, refer to figure 5. 


One Pratt & Whitney J57-P7 Vith afterburner 

15,000 

9,220 

8,000 


19,662 

24,800 


29.5 

29.5 


1 * 
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Figure 1.- Three-view drawing of airplane with original vertical tail and 
the extended as well as the original wing. All dimensions in inches. 
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Figure 2.- Photograph of the airplane. 



E** X 622 

Figure 3 .- Photograph of two airplanes showing tails A and C. 
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Figure 4.- Sketch of vertical tails A, B, and C. Refer to table I for 
physical characteristics of the vertical tail. 
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Figure 5*- Approximated variation of principal moments of inertia and 
inclination of principal axis relative to the body axis. Clean 
configuration. 
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h p , ft 


1.5 
1.4 
33 xIO 3 
32 


a, /3, deg 


p, r, radians/sec 
°t, Q 


4 

2 

0 

-2 

2 

0 

-2 



0 


6 

t, sec 


8 
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(a) Configuration C. M « 1.45; hp *» 52,500 feet. 

Figure 6.- Time histories of lateral oscillations induced by a rudder 

pulse. 
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(b) Configuration D. M 0.78; hp » 40,400 feet. 
Figure 6. - Concluded. 
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Film scale factors per Inch 

p = A»19 radians/sec 
r = 0. 510 radlans/sec 
(3 = IO .75 deg 
a t = 2.50 g 



iPhl _ 
| r bl 


11.7 


x ^ . 19 

10.75 


X 57.3 


7.08 



2.82 

11.7 


2.50 

10.75 


x 57.5 = 2.95 


Figure 7 • - Typical use of the semilogarithmic sheets for determining 
and amplitude ratios, hp = 30,280 feet; M = 0.775; P = 2.98 sec; 
T 1/2 = 2. kO sec. 
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0,02k 

0.068 


Corrected phase angles 

«p b e * \ fSl - V, * - 215 - 1 ' 5 -° = ' 220 - l ° 

»r b |3 = *r b U i - «|) Pl - -75-1 -5.0 = -80.1° 

••HP = *•»!»! - ®PP1 = 2 °5-9 -5-° - 198.9“ 

(a) Determination of phase lag and magnification factor of 0-vane due 

to vane location. 

Figure 8. - Typical sequence employed in the determination of lateral 
derivatives using flight data and the time-vector method. 
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(b) Determination of correct value of Correction necessitated by- 

lateral accelerometer location. 


Figure 8.- Continued. 
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. 2x11! . 2 T Ld t 
2 Id 2 id 


M 

IpI 


+ C, 


L t2 1 


P I 


= 0 



(c) Determination of 


ill 

IPT 



and. positions of the 


tors relative to 3. 


r, 


\|r, and cp vec- 


Figure 8.- Continued. 
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- C Y r 


(°Y r " Cy«) 





Cy p = -0.573 


C-y — C y 
If * i 


. = ( C Yr» ' C Y6 } T1 _ 0.004 _ 

P lr« 1 0.0498 

IP I 


0.080 


(d) Determination of Cyp and ^Cy^ - Cy^. 
Figure 8. - Continued. 
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(e ) Determination of C n p and (C nr - Cnjjj * 


Figure 8.- Continued. 
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Jx ]_£_[ _ Jxg jjJ 

qSb X I p I qSb X I P I 


i£ 

IP 


c l« - c l. 


Ip' I 

FT 


= 0 



c z = -0.0447 



Ip' I 0.151 

IP I 


(f) Determination of and C 7 . 

6 P L V 

Figure 8 .- Concluded. 
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Incompressible dynamic pressure, Ib/sq ft 


Figure 9 . - Effect of flexibility of vertical tail C on vertical tail con- 
tributions to C np , Cn r , C Yp , and C Yr - 6 r = 0 °; a = 0 °; p = -\|r. 

(Manufacturer's estimate.) 
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Figure 10.- Change in lateral stability derivatives due to torsional 
flexibility of the vertical stabilizer in configurations C and D. 
Estimates based on calculated derivatives. 
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(a) hp = 40,000 feet. 


(b) hp = 30,000 feet. 
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Figure 11.. 


(c) h p = 20,000 feet. 

Estimated weight rate of air to jet engine to maintain 
cruising speed. 
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Figure 12.- Estimated contribution of the intake air of the jet engine 

to C Yn and C n . 
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(a) M = 0.80. 


(b) M = 1.20. 


Figure 13.- Influence of C np and 

at M = 0.80 and M = 1.20 at 
curves. 


C l r on "the lateral stability derivatives for configuration C 
h p = 40,000 feet. Solution based on data obtained from faired 
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Condition 

c np 

. Determined Derivatives 

C n 

n P 

c n r ” c np 

1 

O .0985 

- 0.0271 

-0.540 

2 

0.09914- 

- 0.0156 

-0.299 

3 

0.0974 

-0.0432 

- 0.378 



32 Hi . 
qSb X 1 p i ~ 



I r ' 




Determined Derivatives 

Condition 

c n p 

% 

c_ - c„. 

n r n p 

1 

0.0090 

-0.0045 

-0.178 .. 

2 

0.1001 

0.0105 

- 0.148 

3 

0.0980 

-O.OI 55 

-O.I 96 


(b) M = 1.20. 


Figure 14.- Influence of ±1 percent change in C n p on the lateral sta- 
bility derivatives C n ^ and ^Cn r - C n p^ for configuration C at 
M = 0.80 and M = 1.20 at h p = 40,000 feet. Basic values of C n(3 

used were obtained from faired flight data, vector solutions of C^, 
and equation ( 23 ). 
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ij 

V 
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:ermined derivatives 

°ip 

C, 

l P 

% 

C n r -Cnp 

1 

rl 55 « 7 ° 

-0.065 

-0.240 

^0.0980 

- 0 . 34.1 

2 

k60.7 o 

-0.062 

-0.198 

O.O976 

- 0.379 

? 1; 

(-150.7 0 

-O.O69 

-0.290 

0.1000 

-0.283 


] 5 ~} (C n ■ C n* 

*1 J n r n s 


1 r» 

IP 


S * ill - -0.00396 


rlc x rir! = 0.0(4.56 A 


qSb | p | 



Conditior 

a 

pr 

C, 

etermin 

ad^derlvatives 

Clp 

% 

-*r 

5 n r " c np 

1 

-142 0 

-O.O53 

-0.334 

0.0973 

-0.188 

2 

q 

-0.050 

- 0.291 

0.0968 

-0.200 

3 

“ 157 ° 

-O.O58 

- 0.393 

0.0973 

- 0.165 


(a) M = 0.80. 


(b) M = 1.20. 


Figure 15.- Influence of 15° change in <*>p r on the lateral stability derivatives for configura- 
tion C at M = 0.80 and M = 1.20 at h p = V0,000 feet. 
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Condition 

*d 


C n r " C nfi 

1 

AJ.° 

w ■ — 


-0.339 

2 

5. 9° 


-Q.323 

3 

fot 9 s . - 


-0.331 . 


(a) M = 0.80. 



Condition 


*d 


° n r ’ C n* 

1 


3.7° 


-0.183 

2 


3.2° 


-0.170 

3 


! 


-0.200 


(b) M = 1.20. 


Figure l6.- Influence of +0.5° change in on the lateral stability 

derivative ^C nr - C n j^ for configuration C at M = 0.80 and M = 1.20 
at hp = 40,000 feet. 
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Figure 17*- Variation of trim angle of attack with Mach number. 
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(a) Period characteristics. 

Figure 18.- Period and damping characteristics of the airplane as func- 
tions of Mach number, altitude, and angle of attack per load factor. 
Configuration A. 
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(b) Damping characteristics. 


Figure 18.- Concluded. 
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(a) Amplitude ratio characteristics. 

Figure 19.- Amplitude ratio and phase angle characteristics of the air- 
plane at its natural frequency as functions of Mach number and alti- 
tude. Configuration A. 
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(b) Phase angle characteristics. 
Figure 19.- Concluded. 
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(a) Static derivatives. 

Figure 20. - Variation of static and dynamic lateral stability derivatives 

with Mach number. Configuration A. 
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(b) Dynamic derivatives. 
Figure 20.- Concluded. 
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(a) Period characteristics. 

Figure 21.- Period and damping characteristics of the airplane as func- 
tions of Mach number, altitude, and angle of attack per load factor. 
Configuration B. 
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(b) Damping characteristics. 
Figure 21.- Concluded. 
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(a) Amplitude ratio characteristics. 

Figure 22.- Amplitude ratio and phase angle characteristics of the air- 
plane at its natural frequency as functions of Mach number and alti- 
tude. Configuration B. 
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(b) Phase angle characteristics. 
Figure 22.- Concluded. 
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(a) Static derivatives. 

Figure 23.- Variation of static and dynamic lateral stability derivatives 

with Mach number. Configuration B. 
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(b) Dynamic derivatives. 


Figure 23.- Concluded. 
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(a) Period characteristics. 

Figure 2k. - Period and damping characteristics of the airplane as func- 
tions of Mach number, altitude, and angle of attack per load factor. 
Configuration C. 
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(b) Damping characteristics. 
Figure 2k. - Concluded. 
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(a) Amplitude ratio characteristics. 

Figure 25.- Amplitude ratio and phase angle characteristics of the air- 
plane at its natural frequency as functions of Mach number, altitude, 
and angle of attack per load factor. Configuration C. 
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(b) Phase angle characteristics. 


Figure 25.- Concluded. 
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(a) Static derivatives. 

Figure 26.- Variation of static and dynamic lateral stability derivatives 

with Mach number. Configuration C. 
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(b) Dynamic derivatives. 
Figure 26.- Concluded. 
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Figure 27.- Period and damping characteristics of the airplane as func- 
tions of Mach number , altitude , and angle of attack per load factor. 
Configuration D. 
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(b) Damping characteristics. 


Figure 27.- Continued. 
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(c) Damping characteristics. 


Figure 2J.- Concluded. 
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(a) Amplitude ratio characteristics. 

Figure 28.- Amplitude ratio and phase angle characteristics of the air- 
plane at its natural frequency as functions of Mach number, altitude, 
and angle of attack per load factor. Configuration D. 
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(b) Phase angle characteristics. 


Figure 28.- Concluded. 
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(a) Static derivatives. 

Figure 29. - Variation of static and dynamic lateral stability derivatives 

with ifech number. Configuration D. 
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(b) Dynamic derivatives. 
Figure 29.- Concluded. 
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(a) Period characteristics. 

Figure 30.- Summary of period and damping characteristics of configura- 
tions A, B, C, and D as functions of Mach number at h p = h0,000 feet; 
a n = 1.0. 
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(b) Damping characteristics. 
Figure 30. - Concluded. 
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Figure 31.- Summary of the characteristics of the amplitude ratios 
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as functions of Mach number at hp = lj-0,000 feet; a n = 1.0. 
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(a) Static derivatives. 

Figure 32.- Summary of the lateral stability derivative characteristics 
of configurations A, B, C, and D as functions of Mach number at 
hp = 10,000 feetj a n = 1.0. 
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(t>) Dynamic derivatives. 


Figure 32.- Concluded 
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(a) Static derivatives. 

Figure 33.- Comparison of the variation of the lateral stability deriva- 
tives of configurations C and D with Mach number as determined from 
flight with calculated variation, hp = 40,000 feet; a n = 1.0. 
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(b) Dynamic derivatives. 
Figure 33. - Concluded. 
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Figure 34.- Period, damping, amplitude ratio and phase angle 

characteristics as functions of angle of attack at constant Mach num- 
ber. Configurations C and D. 
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Figure 34.- Concluded. 
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Figure 35* - Static and dynamic lateral stability derivative characteris- 
tics as functions of angle of attack at constant Mach number. Con- 
figurations C and D. 
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